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摘要 : ВМА 二 级 结构 比 一 级 结构 包含 更 多 的 信息 , 在 物种 系统 发 育 分 析 中 更 能 反应 真实 情况 。 因 而 ,本 研究 对 得 
HRH Bactrocera correcta 和 瓜 实 蝇 В. cucurbitae 的 L-rRNA 全 基因 进行 了 测序 , Е ГОЙ Н L-rRNA 基因 二 级 
结构 的 模式 图 , 并 分 析 了 其 结构 特征 , 而 后 基于 H45 ~ Н47 荃 环 结构 参 数 和 结构 序列 及 L-rRNA 结构 序列 对 双 怒 
H 13 科 的 系统 发 育 关 系 进 行 了 初步 分 析 。 结 果 表 明 : 双 翅 目 昆 虫 具 有 保守 的 L-rRNA 二 级 结构 ; 不 同 结构 区 碱 基 
组 成 和 分 布 不 均匀 , 区 和 V 区 的 完全 保守 碱 基 含 量 最 高 ，CC 含量 最 高 , 绝 大 部 分 科 YI 区 的 AT 斜率 均 小 于 0; 
科 内 特有 碱 基 和 科 间 少数 保守 碱 基 大 部 分 为 С 或 C。 疼 蚊 科 与 双 怒 目 其 他 科 具 有 较 远 的 系统 发 育 关 系 ， 而 丽 蝇 
科 、 寄 蝇 科 和 食 蚜 蝇 科 有 相近 的 系统 发 育 关 系 ; URS RA Ea ж; 蝶 科 与 蚊 科 同 在 一 个 大 的 系统 
发 育 分 支 。 利 用 单方 面 参数 不 易 得 到 理想 的 进化 结果 , 准确 系统 发 育 分 析 需 要 结合 多 方面 的 参数 。 
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Sequencing of mitochondrial L-rRNA genes of Bactrocera cucurbitae and 
B. correcta ( Diptera: Tephritidae ) and the analysis of structural 


characteristics and phylogenetic relationship of L-rRNA genes in dipterans 
WU Pei-Fu*, XIONG Zhong-Ping*, XU Lei, HUANG Yang-Yang, РАМ Yong-Zhi* ( Key Laboratory of 
Forest Disaster Warning and Control in Yunnan Province, Southwest Forestry University, Kunming 
650224, China) 

Abstract: The RNA secondary structure contains more information than its primary structure, and can 
indicate the phylogenetic relationship more really. Therefore, the complete sequences of mitochondrial L- 
rRNA genes of Bactrocera cucurbitae and B. correcta were sequenced in this study, and the structural 
model of L-rRNA genes in Diptera was constructed. The characteristics of L-rRNA secondary structure 
was analyzed, and preliminary phylogenetic features were inferred based on stem-loop parameters and 
structural sequences of H45 – H47, and structural sequences of L-rRNA genes among 13 families in 
Diptera. The results showed that the dipterans share the conserved secondary structure of L-rRNA gene; 
the nucleotide composition and distribution vary within different structure regions of L-rRNA gene; the 
highest content of completely conserved nucleotides was found in region М and V ; the region М and V 
also contains the highest GC content; the AT-skew is negative in region VI of L-rRNA gene of the 
majority of dipteran families; the unique bases within individual family and individually conserved 
nucleotides between families present higher content of G or C. The family Cecidomyiidae exhibits much 
farther distance with other dipteran families, while Calliphoridae, Tachinidae and Syrphidae share much 
closer phylogenetic relationship; Tabanidae and Nemestrinidae belong to a small branch ; Ceratopogonidae 
and Culicidae share a large phylogenetic branch. The satisfied phylogenetic result would not be obtained 


easily using single phylogenetic parameter, and it is thus necessary to infer accurate phylogenesis through 
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combined parameters of different features. 
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GA Н ( Diptera) &uf& lc, №, 3]. do, ВАЗЕ, 是 
种 类 最 多 的 昆虫 目 之 一 ， 呈 全 球 分 布 , 包括 172 ~ 
179 Ж}, 27 240 000 种 ,其 中 被 报道 的 种 约 为 
132 000#}( Skuhravá and Roques, 2010). ХХ H 
中 , 某 些 昆虫 具有 重要 的 医学 和 社会 学 意义 ,如 ， 
蚊子 [ 蚊 科 ( Culicidae) ] 作 为 载体 可 传播 症 疾 、 登 划 
热 、 西 尼罗河 热 、 黄 热 病 、 脑 炎 等 重大 传染 病 ; 某 
些 昆虫 对 农业 、 林 业 和 畜牧 业 有 重大 影响 , 如 , Ж 
实 蝇 属 Bactrocera 昆虫 可 寄生 于 果蔬,， 对 农业 造成 
巨大 的 经 济 损失 ; 某 些 昆 虫 属 于 益 虫 , 可 寄生 于 或 
捕食 其 他 害虫 , 或 可 传播 花粉 (Roques, 2010) ， 如 ， 
头 蝇 是 绝 大 多 数 寄 生性 害虫 的 天 敌 。 由 此 可 见 , xx 
3H H BBA ACSA BE, H Beak la) Bead 
啊 人 类 生活 , 具有 重要 的 研究 价值 , 是 科学 研究 史 
上 必 不 可 少 的 一 部 分 。 

ЕЖ, ВРЕТ ЛЕ, ХОЗ A Е 
虫 分 子 生物 学 方面 的 研究 得 到 了 迅速 发 展 。 线 粒 体 
DNA 作为 研究 材料 具有 其 独特 的 优越 性 , 成 为 研究 
系统 发 育 的 重要 途径 ( 张 亚 平 和 施 立 明 , 1992) 。 国 
内 外 对 双 翅 目 少数 种 类 的 线粒体 全 基因 组 进行 了 测 
序 , 并 在 基因 组 学 基础 上 做 了 比较 分 析 和 系统 发 育 
人 研究。 目前 , 在 GenBank 中 可 查阅 到 28 FP HGH Н 
昆虫 的 线粒体 全 基因 组 序列 ， 如 黑 腹 果 晶 
Drosophila melanogaster ( Lewis et al., 1995; Ballard, 
2000). IR № fA Aedes aegypti Fil SX РЕ ЖИ Culex 
quinquefasciatus ( Behura et al., 2011). ЖИ 
Bactrocera oleae ( Nardi et al., 2003) #®„ Wiegmann 
等 (2011) 通 过 30 kb 核 基因 、 线 粒 体 全 基因 组 序列 
和 昆虫 形态 学 特征 , 对 149 科 昆 虫 的 遗传 系统 发 育 
关系 进行 了 人 研究 , UNNA LH Bb HUS. 3 次 快速 进 
化 。 此 外 , Ая A AL {ЖЕН 22 0А р FE 
ЭЯ] BA ASR AKA ST FE. 557; DR E 
(2009 ) 通过 线粒体 ND1 和 16S rRNA 基因 序列 探讨 
ТР Ж 12 科 的 系统 发 育 关 系 。 然 而 ,生物 大 
分 子 的 高 级 结构 与 功能 关系 更 加 密切 , 日 其 二 级 结 
构 在 系统 发 育 过 程 中 具有 更 高 的 保守 性 , 更 能 反映 
物种 发 育 系统 真实 的 系统 发 育 史 ( Dai and Wang， 
2008), AX WEE Xt BH A 19 SE WR Bactrocera correcta 
Bezzi 1 А3508 В. cucurbitae Coquillett 的 L-rRNA Ж 
(165 ВМА 基因 ) 进行 了 测序 , 通过 与 GenBank 


中 其 他 双 翅 目 昆 虫 L-rRNA 基因 的 比 对 分 析 , 构建 
J OGG BB He L-rRNA 基因 的 二 级 结构 图 ,并 对 其 
结构 特征 和 系统 发 育 关 系 做 了 分 析 , 从 而 为 双 翅 目 
L-rRNA 基因 结构 系统 发 育 关系 研究 左 定 理论 基础 。 


1 材料 与 万 法 


1.1 ix 

AS Dt Е Br Fi a SE № AJ ЗЕ ЖАТА Г ( 2 10 
3) ARE TS T 2010 年 7 月 采 目 云南 省 元 江 
Яо ЖЕНУ, 将 采集 的 实 蝇 带 回 实 验 
室 , 并 通过 外 形 、 翅 纹 、 翅 脉 、 触 角 等 形态 特征 进 
行 仔细 鉴定 , 将 鉴定 后 的 实验 样品 用 95% 酒精 固 
E, ET -20 冰箱 保存 备用 。 
1.2 基因 组 DNA 提取 与 纯化 

ЖН А ЗЕ Ў Л ЗЕ ШЕ, 利用 基因 组 
DNA 提取 试剂 盒 ( Qiagen, 德国 ) 提取 实 蝇 基 因 组 
DNA。 具 体操 作 步 又 参考 试剂 盒 说 明 书 。 配 置 1% 
琼脂 粮 凝 胶 , 检测 提取 DNA 的 完整 性 。 将 DNA BE 
板 置 于 -207 冰 箱 保 存 备用 。 分 别 用 工头 实 蝇 的 模 
OXY Ж С А BT ЗЕ WY L-rRNA 基因 进行 
测序 。 
1.3 引物 设计 、PCR 扩 增 及 产物 测序 

引用 徐 浪 等 (2007 ) 设 计 的 第 17 ~ 19 对 引物 来 
扩 增 实 蝇 线粒体 L-rRNA ÆA, Вр: N1-J-12314 和 
LR-N-13176; LR-J-12887 和 LR-N-13780; LR-J- 
13417 和 SR-N-14588 。 引 物 由 上 海 生 工 生 物 有 限 公 
司 合成 。PCR 反应 总 体系 为 50 uL, 其 中 2 x Power 
Taq PCR MasterMix (BioTeke ) 25 uL, DNA 模板 
2 uL, 20 nmol/L 上 和 下游 引物 各 1 uL, KK 
21 uL, PCR 扩 增 条 件 为 94% 预 变 性 120 s, H 30 
个 循环 。 每 个 循环 程序 为 : 94°C ЗЕЕ 45 s, 55°C 退 
Х 45 s, 72°С 延伸 1 min 10 s, 最 后 72% 总 延伸 
4 min, 1% ЗАДНЕЕ РСК 扩 增 产物 效果 。 
TI FA SHS BH BE Be DNA 回收 试剂 盒 ( BioTeke ) 回收 
PCR 产物 ,而 后 送 交 北京 百 泰克 生物 技术 有 限 公司 
直接 双 辣 测序 。 
1.4 序列 分 析 

将 双 回 所 测序 列 在 DNAMAN V6.0 中 进行 比 
对 ,结合 测序 波峰 图 确定 所 测序 列 的 质量 。 如 果 测 
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FARD, 将 PCR 产物 进行 重新 测序 。 利 用 
GenBank 中 Blast 程序 来 确定 所 测序 列 的 准确 性 。 
通过 MEGA V5.0 ЖУТО Н РГ L-rRNA 基因 
НЕ SEE ш, 并 计算 各 科 AT 斜率 (AT- 
skew ) 和 完全 保守 核 音 酸 百 分 比 ( 和 INUC ) 的 平均 
E, 其 中 , AT RPE = (A 含量 - ТЎН) /(А 含量 + 
T rt) 。 序 列 比 对 时 所 引用 的 基因 序列 见 表 1, 共 
13 科 38 种 。 

1.5 L-rRNA 二 级 结构 建 模 与 系统 发 育 分 析 

与 以 往 报 道 的 二 级 结构 进行 比 对 , ЖН ЖН ЛИ 
Scie L-rRNA 基因 序列 , 通过 Mfold 程序 (http:// 
www. bioinfo. rpi. edu/applications/mfold/ ) 来 构建 双 
WH L-rRNA 的 二 级 结构 模型 。 模 型 构建 后 , 利用 
Clustal W 程序 比 对 双 翅 目 L-rRNA 序列 ,而 后 在 模 
型 图 中 标注 出 完全 保守 的 核 苷 酸 、 实 蝇 科 内 特有 的 
核 音 酸 、 绝 大 部 分 保守 的 核 苷 酸 及 科 间 个 别 保守 的 
IZER. 

参考 Caetano-Anollés( 2002) 的 方法 , Ж) HAH 
H45 ~ H47 的 7 个 结构 特征 参数 对 双 翅 目 科 间 的 系 
统 发 育 关 系 做 了 分 析 : ® 16 (5), ERA 
(L)、 连 接 桥 中 可 其 代 碱 基数 (B)、 蔡 上 可 符 代 碱 
FERC). BLAZE ERA 98 2829 СТ). BEA ERP 
中 特有 碱 基数 (U) 和 葵 环 大 小 (也 ) 。 对 于 每 个 二 级 
结构 来 说 ,这 些 参数 可 构成 相应 的 集合 Q = 1S, 
L;, B;, Cj, T;, U;, Hj , 其 中 ， 1=1, 2, -,n, A 
S; =5,, L; =L, B; =В,, С, =С,, T; ЕТ,, U;-U, 
Н, =н, 那么 ,可 以 说 两 个 二 级 结构 是 完全 保 
守 的 。 

此 外 , 为 了 进一步 验证 结构 系统 发 育 关系 ,本 
研究 利用 Mega 5.0 软件 分 析 了 基于 L-rRNA 全 序列 
和 H45 ~ H47 结构 序列 的 系统 发 育 关系 ,并 做 了 
比较 。 

1.6 L-rRNA 结构 序列 的 系统 发 育 信号 分 析 

利用 Clustal W 程序 对 双 翅 目 线粒体 L-rRNA Ж 
结构 系列 进行 比 对 分 析 , 并 按照 结构 特征 手动 矫 
正比 对 结果 。 而 后 , 利用 MEGA V5.0 来 分 析 双 起 
目 科 间 的 p ERES. Kimura 双 参 数 距离 、 平 均 转 换 率 
和 平均 颠 换 率 。 计 算 p 距离 和 Kimura 双 参 数 距 离 
HY, Bootstrap 值 为 1 000, 采用 参数 模型 为 “s +v+ 
С”, TARR ЕН], 采用 “GTR +С + Г 
模型 。 


2 fA 


2.1 WWE L-rRNA 基因 特征 

根据 L-rRNA 二 级 结构 特征 , 本 研究 对 各 科 不 
同 分 区 ( 共 6 区 ) 核 车 酸 的 碱 基 售 量 、AT 笠 率 (AT- 
skew) 和 完全 保守 核 车 酸 百 分 比 (%INUC) 进 行 了 分 
析 ( 表 2)。 本 研究 中 ,所 测 番 石 榴 实 蝇 和 瓜 实 蝇 L- 
rRNA 基因 的 大 小 分 别 为 1 326 bp 和 1 327 bp. 
2.1.1 HH L-rRNA 基因 各 结构 区 GC 含量 : 由 
d 2 ИЖ, ЙН L-rRNA 基因 各 结构 区 的 碱 
基 镶 量具 有 差异 性 。 从 整体 上 来 看 ，G + C% 的 大 
小 顺序 为 No» Vo» M. fH EE SCR 
( Cecidomyiidae) , RIF HAP} ( Syrphidae ) ЖП] Р 
( Nemestrinidae) МХ G + C% & E Pob E IX, 
HAP T KH С + C96 含量 (30. 94% ) BERT 
Hd 44 Ko Xp SE № Ph ( Tephritidae ). № F} 
( Culicidae ), № fl ( Ceratopogonidae ), tÈ f 
( Tabanidae ). Е № fL ( Oestridae ) ЖОЮЫ} 
( Drosophilidae ) , # fl ( Muscidae ) 1E WB 
( Agromyzidae ) ЧР В HEB} ( Syrphidae) № I, ПЖ 
亚 区 来 说 , 其 G+ Coo 含量 的 大 小 顺序 为 : I > П 
> Ш. #7 № $} ( Tachinidae ) Ж py x8 mr B 
( Nemestrinidae) ) G +C% Ge Ш> I > Il; my ee 
fi ( Calliphoridae ) МХТ: ( Cecidomyiidae ) С + 
С% EE I» Ш> П. 58 H № L-rRNA 基因 各 
结构 区 AT 含量 变化 与 GC 含量 刚好 相反 , 这 里 不 
再 做 说 明 。 
2.1.2 HE L-rRNA 基因 各 结构 区 AT 斜率 : 在 
ХОЙ SPP, АТ 碱 基 分 布 特征 也 具有 明显 差异 
性 。 整 体 而 言 ， 除 狂 蝇 科 ( Oestridae ) ЖП 3H at fl 
( Метеѕітіпідае) К, МІХ AT 斜率 均 小 于 0。 不 同 科 
间 AT PFERD ЖАН, МНР, БР. 
А. ВРТА ТАВЬ У 区 的 АТ 斜率 高 于 其 
他 各 结构 区 ; ВВ, НОР, ap RB A ЖШ БЕШ 
区 的 AT 斜率 高 于 其 他 各 结构 区 ; BOR. В Ж 
РЕР СВ I 区 的 АТ 斜率 高 于 其 他 各 结构 区 ; S 
BE I KKJ АТ 斜率 高 于 其 他 各 结构 区 。 特 别 指出 
MHP I, Ш, VA VL PERS AT 斜率 均 为 负 值 。 
从 各 科 整 体 AT 斜率 来 看 , 所 有 和 斜率 均 大 于 0。 
2.1.3 SGA L-rRNA 基因 各 结构 区 % INUC 的 特 
ÎE: 由 表 2 可 以 看 出 , ЙН L-rRNA 基因 中 , МК 
和 V 区 的 %INUC 值 最 高 , 均 大 于 其 他 各 结构 区 的 
统计 值 ; 除 果 蝇 科 工区 和 亚 区 的 统计 值 相同 外 , 其 
他 两 科 亚 区 的 %INUC 值 最 小 。 从 科 间 角度 来 看 ， 
条 蝇 科 的 %INUC 值 最 大 , 实 蝇 科 与 蚊 科 的 %INUC 
值 近 似 。 


924 


科 
Family 


Tephritidae 


Culicidae 


Calliphoridae 


Ceratopogonidae 
Tabanidae 


Oestridae 


Drosophilidae 


Tachinidae 
Muscidae 


Agromyzidae 


Cecidomyiidae 


Syrphidae 


Nemestrinidae 


Table 1 
属 


Genus 


Bactrocera 


Ceratitis 


Culex 


Aedes 


Anopheles 


Chrysomya 
Cochliomyia 
Lucilia 
Culicoides 
Cydistomyia 
Dermatobia 
Hypoderma 


Drosophila 


Exorista 
Haematobia 


Liriomyza 


Mayetiola 
Rhopalomyia 


Simosyrphus 


Trichophthalma 
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表 1 系统 发 育 分 析 中 所 涉及 的 序列 


The sequences used in the phylogenetic analysis 


种 
Species 
Bactrocera cucurbitae 
B. carambolae 
B. correcta 
B. dorsalis 
B. minax 
B. oleae 
B. papayae 
B. philippinensis 
B. tryoni 
Ceratitis capitata 
Culex quinquefasciatus 
C. pipiens pipiens 
Aedes aegypti 
A. albopictus 
Anopheles darlingi 
A. gambiae 
A. quadrimaculatus A 
Chrysomya putoria 
Cochliomyia hominivorax 
Lucilia sericata 
Culicoides arakawae 
Cydistomyia duplonotata 
Dermatobia hominis 
Hypoderma lineatum 
Drosophila littoralis 
D. mauritiana 
D. melanogaster 
D. sechellia 
D. simulans 
D. yakuba 
Exorista sorbillans 
Haematobia irritans 
Liriomyza sativae 
L. trifolii 
Mayetiola destructor 
Rhopalomyia pomum 
Simosyrphus grandicornis 


Trichophthalma punctata 


GenBank 登录 号 


GenBank accession no. 


NC_016056 
NC_009772 
JX_456552 
NC_008748 
NC_014402 
NC_005333 
NC_009770 
NC_009771 
NC_014611 
NC_000857 
NC_014574 
NC_015079 
NC_010241 
NC_006817 
NC_014275 
NC_002084 
NC_000875 
NC_002697 
NC_002660 
NC_009733 
NC_009809 
NC_008756 
NC_006378 
NC_013932 
NC_011596 
NC_005779 
NC_001709 
NC_005780 
NC_005781 
NC_001322 
NC_014704 
NC_007102 
NC_015926 
NC_014283 
NC_013066 
NC_013063 
NC_008754 





NC_008755 


序列 长 度 
Size (bp) 
1 327 
1 324 
1 326 
1 333 
1 333 
1 323 
1 325 
1 325 
1 326 
1 335 
1 333 
1 333 
1 335 
1 333 
1 394 
1 325 
1 321 
1 329 
1 324 
1 327 
1 314 
1 337 
1 324 
1 316 
1 324 
1 323 
1 325 
1 322 
1 324 
1 326 
1 329 
1 321 
1 323 
1 321 
1 302 
1 264 
1 339 
1 318 
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2.2 MWH L-rRNA 基因 二 级 结构 模型 构建 与 系 
统 发 育 特 征 

利用 Lasergene 中 MegAlign 程序 分 析 表 明 ， 双 
3H H rP3E 453 个 完全 保守 核 车 酸 , HH, A н 
33.77% , T 41. 5096 , С 15.45%, C 9. 2795, 
完全 保守 碱 基 中 ，G 含量 显著 高 于 双 翅 目 各 科 L- 
rRNA 基因 的 平均 С 含量 ( 表 2) 。 结 合 表 2 AA 1, 
可 以 看 出 绝 大 部 分 完全 保守 核 车 酸 位 于 VY 区 和 V 
区 。 大 部 分 荃 环 结构 上 的 碱 基 在 各 科 内 呈 高 度 保 
SP. 但 这 些 碱 基 在 科 间 具有 较 高 的 颠 换 率 ( 即 A 和 
T HEH). 
2.2.1 ХН L-RNA 基因 二 级 结构 模型 . 除 达 氏 
按 蚊 Anopheles darlingi L-rRNA 基因 的 大 小 为 1 394 
bp 外 , 其 他 双 翅 目 昆虫 的 L-rRNA 基因 大 小 均 约 为 
1 326 bp, 且 具 有 81% 以 上 的 相似 度 , 这 为 双 翅 目 


L-rRNA 基因 二 级 结构 模型 的 构建 提供 了 便利 。 双 
ЙН L-rRNA 具有 非常 保守 的 6 个 结构 区 (I ~ 
М), 3551 个 葵 环 结构 , 其 中 ， 工 区 包含 8 FAH 
结构 (图 1) 。 跟 其 他 节 胶 动物 一 样 , 双 翅 目 昆 虫 L- 
rRNA 基因 亚 区 缺乏 明显 的 荃 环 结构 。 从 图 中 可 以 
看 出 ,大 部 分 结构 碱 基 都 呈 高 度 保守 ， 部 分 碱 基 呈 
完全 保守 状态 , 且 大 部 分 保守 碱 基 在 结构 中 均匀 分 
布 。 双 翅 目 部 分 葵 环 结构 呈 高 度 保守 : YR Н23, 2E 
环 H28. № 34. Ж 3$、 环 40、 环 H44 MAH HAT, 
ЖЕНЕР HAT 在 双 翅 目 中 完全 保守 。 从 双 翅 目 科 
内 特有 碱 基 和 科 间 个 别 保守 碱 基 来 看 , 这 些 碱 基 具 
有 整体 特征 ， 即 绝 大 部 分 碱 基 为 C 或 C, Н С 碱 基 
占 绝 大 部 分 ( 如 实 蝇 科 和 果 晶 科 )。 然 而 也 有 例外 ， 
ЗАВ ( Cecidomyiidae) 2H H45 ~ H47 中 特有 碱 基 
HART, 且 A 和 T 的 比例 大 致 相等 。 






т 
-$ ets (09000, 
Н H328- 0% 9 T 
1 os- Ф А ^ Pd 
нове 6% D ooe? 
99 mi $4 H33 2 
тө e 00 oo? % H30 909 H34 e 
97000060006 60 009% e ме 0000000°% $ - 
9999 ^ 00 ei 0000? 9e ee 000 ? Ба H39 
H22 $- $e 8 90679900 09%, 
t -ФН27 roe" H35 "98 $- ноб 2 
ү: TT TA 2 s H41 
B dn et H ELO t 
" e E -00 x tite 1 
H wai: НИНЫ 
2 $-3 H26 st ° 
56 2°% mI e 
H QUO OH I 3 ч ü ww 
0 т $ Ны Ш oe 2 2 e 9-9 H43 
e 0 ð- 0090020000000, 009- - $-9 6.8 
а @® өф Н25 09-0 
$ o Hi! 990° 00006-2 00009900, e H36 ee” 0 
» © H17 ,0-000000000 e? eee Ф% 0-0 Ve ee 
өө rij 900990 % H37 2 a 
62000, Sos 1 oseoeseoeeeoeee IV eeeo0o0eoeoeeeoeoee9 ES е, 
nm S й e" 0000000000 © 0908000000 90 00000000000000000" f 
e e e өө ' т 006000660000000009.0 
“8 9 HIS 009.4 $ A 8 e? ee 299- 2000? 
#% Не т ‚Мы arm SH 
: A ATT pe ҮТҮ, 9 = у 
oo ы H Г ое? $, 09008 000000000,0 e "S $ 
e e 9.9 6- өө @- H49 H48 $3 
$00" "ee 8-8 112 H51 ®-Ф " $-9 H45 
e e т" oo € 
900 e- ^^ m | == e? : 
6090900 нм ce во 6-6 
: cmi: ji Fr 
e TA Ф H46 
n $e 500000000 ё A " 
Sede. ө Soa PA at: 
H15 Ф-өө-ө°®®е„°®о„ © 900° e 
of. H 90700, H H47 
NERIS 906 ө e$ 
00 06 *e e о 
H7 өө 909 





图 1 ЖЕН L-rRNA 基因 序列 所 构建 的 二 级 结构 模式 图 
Fig. 1 The secondary structure model of Diptera based on the L-rRNA sequence of Bactrocera correcta 
BUS E REREH Y Sé PRS ER RAFE), SCHR ARR BUCH BR (ТЛВ), AEB ORE BR RETE) 及 科 间 个 
FRE A CE BR HIRAF); 图 2(A) 同 。The special nucleotides are presented as follows: completely conserved nucleotides (red letters with 
blue background), the unique bases of Tephritidae ( white letters with red background ), most conserved nucleotides ( black letters with green 


background) , individually conserved bases among families ( white letters with magenta background). The same for Fig. 2( A). 
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2.2.2 二 级 结构 系统 发 育 分 析 : 为 了 分 析 双 翅 目 
L-rRNA 基因 二 级 结构 的 系统 发 育 特征 , 本 研究 利 
用 Clustal W 程序 对 L-rRNA 基因 的 序列 进行 了 比 对 
MATRE, 并 结合 结构 比较 选取 了 H45 ~ H47 XE 
环 结 构 作 为 二 级 结构 系统 发 育 分 析 的 代表 结构 (图 
2: A) 。 在 构建 参照 模型 结构 时 ， 对 双 翅 目 L-rRNA 
基因 H45 ~ H47 葵 环 结构 的 序列 做 了 模式 分 析 ( 图 
2: А 中 的 Model; 图 2: B) ， 而 后 根据 模式 序列 来 
构建 参照 模型 结构 。 此 外 , ИМЕ ГӘН Н H45 ~ 
H47 茎 环 结构 区 的 二 级 结构 图 ,作为 外 围 结构 参 
数 。 图 中 双 箭 头 所 示 碱 基 为 科 内 可 符 代 碱 基 。 

由 图 2 可 以 看 出 , Н H45 ~ H47 52 


A 


55 Ж 


区 的 主要 不 同 点 体现 在 以 下 几 个 方面 ( 按 模式 图 ) : 
蕉 长 (S2)、 连 接 桥 大 小 (L)、 连 接 桥 中 可 百代 碱 基 
数 (B)、 葵 上 可 蔡 代 碱 基数 (C2 和 C3) , НИЕ 
特有 碱 基数 (T)、 科 内 连接 桥 中 特有 碱 基数 (U) 和 
环 大 小 (H2) 。 由 图 看 出 瘦 蚊 科 的 结构 参数 变化 非 
常 明显 。 通 过 结构 参数 所 做 的 系统 发 育 图 表明 (图 
3: A) ， 鳞 翅 目 作为 外 围 群 单独 位 于 一 个 大 分 文 
INE SIGE BARRA RIL RRR RA, 
单独 位 于 一 个 分 支 ; SERB SATB ob all ERAT 
支 ; FRAP ASR AU RR ORA. 
FHP ЛЛ x МАР. TRR AR BL ЕЛ 
的 系统 发 育 关 系 , 位 于 一 个 小 分 支 。 
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DOA Н L-rRNA 基因 H45 ~ H47 & АНДАХ (А) 与 结构 模式 序列 分 析 (B) 


Fig. 2 Individual stem-loops of H45 - Н47 (А) and analysis of structure motif sequence (В) of L-rRNA genes in Diptera 
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3 WHA L-rRNA 基因 结构 序列 的 系统 发 育 信和 号 
Table 3 The phylogenic signals of the structure sequences of L-rRNA genes in Diptera 


1 2 3 4 5 6 

1 0.12 0.10 0.12 0.12 0.09 
2 0.10 0.03 0.14 0.15 0.05 
3 0.09 0.03 0.10 0.12 0.02 
4 0.10 0.12 0.09 0.18 0.11 
5 0.10 0.13 0.10 0.14 0.11 
6 0.08 0.05 0.02 0.09 0. 09 

7 0. 10 0.07 0. 03 0. 10 0. 07 0. 05 
8 0.09 0.03 0.00 0.09 0.10 0.02 
9 0.10 0.05 0.06 0. 12 0. 12 0.06 
10 0.06 0.09 0.06 0.09 0. 09 0. 05 
11 0.22 0. 20 0.20 0. 24 0. 23 0.20 
12 0. 11 0.09 0.07 0. 10 0. 10 0.07 
13 0. 13 0.09 0. 08 0. 14 0. 12 0.09 
1 0.04 0.04 0.08 0.04 0.04 
2 0.07 0.01 0.09 0.06 0.01 
3 0.06 0.03 0.08 0.05 0.01 
4 0. 04 0. 05 0. 02 0. 12 0.09 
5 0.08 0.09 0.06 0.06 0.06 
6 0.05 0. 04 0.01 0.02 0.05 

7 0.06 0.03 0.00 0. 02 0. 06 0.01 
8 0. 06 0. 03 0.00 0. 02 0. 06 0.01 
9 0. 08 0. 04 0.06 0.06 0.08 0.06 
10 0. 04 0.08 0.05 0.02 0.06 0. 04 
11 0.23 0.18 0.20 0.25 0.23 0.20 
12 0.08 0.06 0. 04 0.06 0.05 0. 04 
13 0.08 0.03 0.01 0. 04 0. 08 0. 02 


> > > sees 


7 8 9 10 11 12 13 
. 12 0. 10 0.12 0.07 0.32 0. 13 0.16 
.07 0. 03 0. 06 0.10 0.27 0.11 0.11 
04 0.00 0. 06 0. 07 0. 27 0. 07 0. 09 
12 0. 10 0. 15 0. 11 0. 36 0. 12 0. 19 
07 0. 12 0. 15 0. 11 0. 34 0. 12 0. 15 
. 05 0. 02 0. 07 0. 05 0. 27 0. 08 0. 11 
0. 04 0. 10 0. 08 0. 32 0. 09 0. 08 
0. 03 0. 06 0. 07 0. 27 0. 07 0. 09 
0.09 0.06 0. 10 0. 23 0. 12 0. 15 
0.07 0.06 0.09 0. 30 0. 09 0. 12 
0.22 0.20 0. 18 0.21 0. 30 0. 39 
0. 08 0.07 0. 10 0. 08 0.21 0. 15 
0.07 0.08 0. 12 0.11 0.26 0. 12 
0.05 0.04 0.04 0. 03 0.09 0.05 0.09 
0.05 0.01 0.02 0.02 0.09 0.04 0.08 
0.04 0.00 0.00 0.02 0.07 0.04 0.08 
0.10 0.08 0.09 0.09 0.11 0.05 0.15 
0.01 0.05 0.07 0.05 0.11 0.07 0.07 
0.04 0.01 0.00 0.02 0.07 0.04 0.09 
0.04 0.04 0.03 0.11 0.05 0.06 
0.00 0.00 0.02 0.07 0.04 0.08 
0.06 0.06 0.01 0.05 0.04 0.07 
0.05 0.05 0.09 0.05 0.03 0.06 
0.20 0.20 0.18 0.25 0.07 0.16 
0.04 0.04 0.08 0.06 0.23 0.10 


0.01 0.01 0.08 0.06 0.23 0.05 


1; BCR} Culicidae; 2: SHER] Tephritidae; 3: МНЯ Calliphoridae; 4; № Ceratopogonidae; 5; Я Tabanidae; 6; 11| Oestridae; 7: 
Я Drosophilidae; 8: ЗЕ Tachinidae; 9: WP} Muscidae; 10: ЖЕ Agromyzidae; 11; Ж}. Cecidomyiidae; 12: В HAA 








Syrphidae; 13; ЖЖ Nemestrinidae. 表 中 正体 数据 表示 双 翅 目 科 间 的 p 系统 发 育 距 离 ; 粗 正 体 数据 表示 科 间 的 Kimura 双 参 数 系 
统 发 育 忠 离 ; 斜体 数据 表示 科 间 核 苷 酸 的 平均 颠 换 率 ; 粗 斜 体 数据 表示 科 间 核 苷 酸 的 平均 转换 率 。The data with standardized form are 


the p-distance between families of Diptera; those with bold standardized form represent the distance from Kimura 2-parameters model; those with 


italic type are the average of nucleotides transversion rate between families; and those with bold italic type are the average of nucleotides 


transition rate between families. 


征 可 能 L-rRNA 基因 的 功能 相关 。 以 往 人 研究 表明 ， 
这 些 结构 区 是 转 肽 酶 的 结合 中 心 , SOS 核糖 体 亚 基 
在 和 绰 晶 合成 过 程 中 通过 肽 出 口 通道 (peptide exit 
tunnel, PET) 与 L-rRNA 的 5 个 结构 域 结合 (Frank et 
al., 1995; Nissen et al., 2000)。 由 此 可 见 , 动物 L- 
rRNA 结构 高 度 保守 的 原因 与 其 功能 有 密切 关系 。 
REA Н Р L-rRNA 基因 的 绝 大 部 分 结构 序列 


高 度 保守 , 但 科 间 或 科 内 亦 存在 明显 的 转换 和 菇 换 
显现 , 说明 不 同 科 间 具 有 各 自 的 突变 和 系统 发 育 历 
程 。 为 了 维持 二 级 结构 和 功能 的 稳定 ， 双 翅 目 L- 
rRNA 基因 的 系统 发 育 受 结构 和 功能 的 双重 选择 压 
Л, 也 就 是 说 ,虽然 存在 核 苷 酸 的 替代 现象 , 但 双 
WH L-rRNA 基因 的 二 级 结构 处 于 高 度 保守 状态 ， 
有 目 与 其 他 市 胶 动 物 或 峭 椎 动物 的 二 级 结构 具有 很 高 
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的 相似 度 ( Cannone et al., 2002; Negrisolo et al., 
2011 ) 。 

有 意思 的 是 , 双 翅 目 各 科 L-rRNA 基因 不 同 结 
构 区 的 碱 基 组 成 和 分 布 具 有 差异 性 , 且 科 内 符 有 碱 
基 和 科 间 保守 的 少数 碱 基 大 部 分 为 C 或 C。Li 等 
(2008 ) XFER 93 种 鱼 的 L-rRNA 结构 序列 进行 了 
分 析 , 认为 L-rRNA 结构 序列 具有 转换 -三 换 偏 问 
Е, 茎 区 核 车 酸 蔡 代 率 低 于 环 区 ; КАННА 
表现 为 A >T, C>G; 茎 区 碱 基 配对 受到 选择 压力 ， 
Жу С> С, Т> А, Ramirez 和 Ramírez( 2010) 对 
蜗牛 L-rRNA ÈA NV KA V 区 进行 了 测序 分 析 ， 结 
果 表 明 大 部 分 芸 区 呈 保 守 状 态 , 而 环 区 易 突 变 。 与 
本 研究 相 比 , 这 些 结果 可 从 本 研究 图 1、 图 2 和 序 
列 分 析 得 到 反应 ， 即 双 翅 目 L-rRNA 结构 序列 在 整 
体 上 也 有 类 似 的 现象 。 然 而 , 这 些 研究 均 未 对 科 内 
或 科 间 少数 保守 碱 基 进 行 分 析 。 双 起 目 科 内 L- 
rRNA 基因 特有 碱 基 特点 能 否 说 明 在 系统 发 育 历 程 
中 该 基因 CRC 含量 有 增加 的 趋势 , RRA THE 
一 步 人 研究 。 昆 虫 线粒体 DNA 的 一 处 明显 特征 是 АТ 
含量 非常 高 , 且 大 部 分 昆虫 的 A ZEA Г T. 
至 今 为 止 , 双 翅 目 线粒体 全 基因 组 的 AT% 范围 值 
H 67. 3% ~ 85. 296 , 平均 值 为 77.1%, 日 AT 斜率 
均 为 正 值 , 其 中 21 ARP (21/37) 的 AT 斜率 低 于 
平均 值 0.036( 本 人 研究 计算 数据 ) 。 此 外 ,结合 完全 
保守 碱 基 中 高 CC 含量 特征 ,可 以 认为 科 内 或 科 间 
少数 保守 碱 基 组 成 特征 可 能 是 核酸 组 成 选择 压力 的 
2096. XAH L-rRNA 基因 VI 区 来 说 ,大 部 分 科 
的 A 含量 低 于 了 T 含 量 , 且 变 异性 最 大 , 其 大 小 与 L- 
rRNA 基因 大 小 有 直接 联系 (如 达 氏 按 蚁 ), 说 明 VI 
区 可 能 是 L-rRNA 基因 与 转 肽 酶 结合 中 心 的 调控 区 
R, 在 系统 发 育 中 其 碱 基 组 成 具有 选择 性 。 特 别 指 
出 的 是 , BERL 工区 的 GC 含量 (30. 94% ) 显著 高 
于 其 他 各 科 , 而 整体 CC 含量 却 低 于 其 他 各 科 , 且 
其 L-rRNA Ж КАМЕ А ЈУ, ВОР 
具有 与 其 他 科 不 同 的 系统 发 育 历程 , 在 系统 发 育 过 
erp SU Ay Sc, 这 点 可 从 本 研究 的 系统 发 育 分 析 结 
AR SHE жї дд» ХР БОИ, 其 L-rRNA 结构 
FRU BR БЕЛП V 区 外 ,AT 斜率 均 大 于 其 他 各 科 的 
对 应 结构 区 , 整体 斜率 大 于 其 他 各 科 , Ау IL RA Г 
科 可 能 也 具有 与 其 他 科 不 同 的 系统 发 育 历程 。 对 疼 
AP A PY aT BL SCN ASR SRL, 还 需 各 
方面 的 考证 分 析 , 包括 形态 特征 、 生 活 史 等 。 
3.2 ” 双 翅 目 昆 虫 科 间 的 系统 发 育 关系 

近年 来 ,人们 采用 不 同 手 段 对 双 翅 目 昆 虫 的 系 


统 发 育 关 系 做 了 大 量 研究 ( 池 军 等 , 2010; Gibson et 
al., 2010) , Han 等 (2002 ) 及 Han 和 Ro( 2005) 的 系 
统 发 育 分 析 表 明 : MEPR REL. РВ 
具有 较 近 的 系统 发 育 关系 , 其 中 丽 晶 科 和 蝇 科 同属 
一 个 分 文 ; SEWER ИАР ( Pyrgotidae ) 位 于 一 个 
系统 发 育 分 文 ; 食 蚜 蝇 科 与 尖 尾 蝇 科 ( Lonchaeidae) 
位 于 一 个 系统 发 育 分 文 。Wiegmann 等 (2011 ) 研究 
表明 ， 丽 蝇 科 和 寄 蝇 科 同 属 一 个 小 分 文 , 与 狂 蝇 科 
和 蝇 科 有 较 近 的 系统 发 育 关 系 , БЕНИН 
同 在 一 个 大 的 系统 发 育 分 文 ; ЗЕ ВБ ТАВ 
(Pyrgotidae) 位 于 一 个 系统 发 育 分 文 ; UCB A ЯЙ 
VRBE INA; 嵌 科 与 蚊 科 同 在 一 个 大 的 系 
统 发 育 分 文 , 即 长 角 亚 目 分 文 。 在 本 研究 中 , 通过 
结构 参数 进行 的 系统 发 育 结果 表明 瘦 蚊 科 与 双 翅 目 
其 他 科 具 有 较 远 的 系统 发 育 关 系 ,， 而 丽 蝇 科 、 寄 蝇 
科 和 食 蚜 蝇 科 有 相近 的 系统 发 育 关 系 , 这 与 以 往 的 
结果 相 一 致 , fH SCORE, EISE. РТ PE В| 
的 系统 发 育 分 析 比 较 混 乱 , 这 可 能 与 结构 参数 的 取 
舍 、 结 构 区 的 选择 、 系 统 发 育 分 析 方 法 等 有 关 。 基 
于 核 苷 酸 序 列 和 H45 ~ H47 结构 序列 系统 发 育 信 和 号 
的 分 析 , 基本 上 取得 了 与 以 往 相 类 似 的 结果 。 尺 管 
如 此 , 单独 利用 L-rRNA 基因 亦 不 能 得 到 理想 的 系 
统 发 育 分 析 结 果 , 因而 , 本 研究 认为 准确 系统 发 育 
分 析 需 要 结合 多 方面 的 参数 , 如 : 多 基因 序列 或 基 
因 组 序列 、 形 态 特 征 、 基 因 结 构 特 征 、 生 活 史 等 。 
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